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ABSTRACT 


A summary of Bragg deflection theory and various approaches 
to direct detection acoustooptic spectrum analysis (AOSA) 
modeling is presented. A suitable model is chosen and 
extended to include the effects of diffraction efficiency, 
transducer efficiency, irradiance profiles of Peet laser 
illumination, aperture size of the Bragg cell, and the acoustic 
attenuation experienced by the acoustic wavetrain generated by 
the input r-f signal. A Fortran program is developed to model 
the AOSA and predict the output image plane intensity profiles. 
A second version of the program includes a time variable per- 


mitting dynamic simulation of the system response. 
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i. LN ese DUCTION 


A. OVERVIEW 

Assessment of emergent technologies leads to a continuing 
focus on acoustooptical techniques for spectral analysis and 
filtering/excision [Ref. 1]. Integrated optics with distinc- 
tive advantages in size and weight easily meet the require- 
ments for radar warning receivers. Bulk wave devices offer 
superior dynamic range, bandwidth, and resolution over inte- 
grated optical devices. Foreseeable applications include 
surveillance and electronic intelligence. 

An acoustooptical spectrum analyzer transforms a wideband 
time domain Signal into an optical spatially detected fre- 
quency domain signal. Removal of narrowband interference 
may be accomplished by optical excision techniques, to selec- 
tively notch filter unwanted spectral components. U. S. 
military applications in this area could greatly benefit 
from the speed of response and high quality interference re- 
jection capabilities possible with these techniques. The 
potential of these techniques was recognized and led to the 
establishment of a joint NAVELEX 350/DARPA program to examine 
these areas. 

The technical objectives of this effort were: 

1. Survey potential applications and analyze technical 

requirements imposed by each. 


2. Analyze alternative optical excision configurations 


and components. 
laph 





3. Develop and procure major system components. 
4, Fabricate and measure alternative broadband excisor 
configurations to verify analytical results. 

5. Fabricate and field test a brassboard model. 

6. Design a preliminary advanced development model. 
(me principal civilian contractor for this project is PROBE 
Systems, Inc. PROBE, assisted by the developmental efforts 
at the Naval Research Laboratory (NRL) and Naval Postgraduate 
School (NPS), was tasked with the attainment of these objec- 
tives. NRL's specific program assistance was the development 
of a photodichroic optical clipper to be evaluated by PROBE. 
The research summarized in this thesis is a continuation of 
research initiated as part of the NPS contribution to the 


development program. 


Bee bHeolS OBJECTIVE 

The objective of this work effort was the attainment of 
an understanding of the principles of direct detection acousto- 
optical spectrum analysis. Previous work at the Naval 
Postgraduate School resulted in the construction by Regan 
[Ref. 2] of a working laboratory acoustooptical spectrum 
analyzer and the subsequent modification of the existing 
equipment by Smith [Ref. 3] to function as an acoustooptical 
Spectrum excisor (AOSE). The acoustooptical spectral excisor 
utilized optical excision, the effect of an object placed 


in that part of the optical path where the undesirable 
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frequency components have been spatially spread. The optical 
excision research included propagation considerations and 
diffraction by the obstacle. The effectiveness of the optical 
excisor in removing interference is dependent upon these 
factors and more fundamentally upon the ability of the de- 
flected beam irradiance pattern to accurately represent the 


input r-f signal spectral components. 


C. ACOUSTOOPTICAL SPECTRAL ANALYSIS 

An acoustooptical spectrum analyzer (AOSA) consists of 
an acoustooptical light deflector (a Bragg cell), a coherent 
light source (a laser), and a detector. An input r-f signal 
excites a piezo-electric transducer mounted on the acoustic 
medium on the Bragg cell. This launches an acoustic wave- 
train the length of the cell which modulates the refractive 
index of the interaction medium by the elasto-acoustic 
process. The modulated interaction material functions as 
a moving diffraction grating for the incident laser light. 
A one-dimensional angular pattern of diffracted beams is 
produced and focused on a photodetector array. The pattern 
replicates the spectral content of the original r-f signal. 
The intensity of the individual beams is directly proportional 
to the r-£ signal component intensity. The deflection 
angles of the individual beams are proportional to the r-f 
Signal frequency component. The detector array photosites 


Video detect the incident beam patterns and, in a square law 
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fashion, convert the spectral information for electronic 
transmission to a display oscilloscope. Figure 1 shows the 
NPS laboratory AOSA. Figure 2 shows a typical oscilloscope 
display output from the system. Korpel [Ref. 4] includes a 
basics discussion on the theory of Bragg cell applications 
for acoustooptical spectrum analysis in an overview on 


acoustooptics. 


D. SPECIFIC GOALS 

There were three initial goals for this work. The first 
was the study of the factors affecting the Bragg cell 
deflection of incident laser light. The second was the 
analysis of the effect of these factors on the replication 
of the input r-f signal spectrum in the deflected beam 
pattern. The third was the development of a Fortran program 
that would model the acoustooptic spectrum analysis process. 

Mm Bragg Diffraction 

The study and analysis of the factors affecting 

Bragg deflection was focused towards an understanding of the 
ability of the deflected beam pattern to faithfully replicate 
the spectral components of the input r-f signal. This 
knowledge was considered essential in order to construct a 
Viable computer simulation model. A study of Bragg diffrac- 
tion and the various approaches to modeling of the deflected 
beam irradiance profiles was undertaken. Previous works at 


NPS by Regan [Ref. 2] and Smith [Ref. 3] were used as 
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Figure 2. Oscilloscope Display of a 70 MHz Monotonic Input 
Signal to the NPS Acoustooptic Spectrum Analyzer 
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Searting points. Bragg diffraction is traditionally theo- 
rized by use of complex coupled mode wave equations. The 
complexity of these equations rendered them unsuitable for 
the first order simulation sought. A study of other 
approaches was conducted to arrive at a method that would 
provide both adequate accuracy and a reasonable level of 
complexity. 

2. Computer Model 

An approach to the modeling of the AOSA was selected 
based on the work of Hecht [Ref. 5]. He shows the image 
plane intensity as the squared magnitude of the convolution 
of the spatial Fourier transform of an illumination function 
and the transform of the converted acoustic wave. Some 
modifications in the form of transmission factors were applied 
by algebraic multiplication. Hecht's illumination factor is 
a spatial weighting function composed of the optical ampli- 
tude profile, the optical aperture size of the Bragg cell, 
and the acoustic attenuation of the interaction material. 

The added transmission factors were the diffraction efficiency 
and the transducer efficiency. 

A Fortran program was developed to model the AOSA 
given the specific physical parameters of the system and 
desired Simulation requirements such as laser wavelength and 
laser beam profile. This program was used to study the 
quality of the replication of the r-f signal spectral components. 


The major thrust of this effort was to produce this model for 


Le) 





utilization with previously developed optical excision 
programs. While the program is configured to simulate 
various component deflectors, 1t was structured about the 
AOD-70 model Bragg deflector, manufactured by the 
IntraAction Corp., which is scheduled for future installa- 
tion in the NPS laboratory AOSA. The AOD-70 Bragg deflector 
is constructed with high quality flint glass and Indium 
bonded Lithium Niobate stepped array transducers. It has 
a 50 MHz to 90 MHz bandwidth and a center frequency of 70 MHz. 
During the development of the program it became 
apparent that the addition of a third dimension, time, to 
the model would significantly enhance the ability of the 
model to demonstrate the response of the AOSA under various 
Simulation conditions. A second Fortran program was written, 
essentially a variation of the first, to provide this dynamic 


capability. 


m. PERSPECTIVE 

The application of both AOSA and AOSE techniques in wide- 
band, high density signal environments is dependent on the 
accuracy of the diffracted beam profile in its replication 
of the input signal spectra. The classical model for 
determination of the diffracted beam profile is presented 
by Hecht [Ref. 5]. In a simplified form the diffracted beam 
profile is the convolution of the spatial Fourier transform 


of the illumination field incident on the Bragg cell and the 





Spatial transform of the generated acoustic wave within the 
cell. In detail both the illumination function and the 
acoustic wave are subject to modifying factors. The signifi- 
cant factors include diffraction efficiency, transducer 
efficiency, the irradiance profile of the incident laser 
light, the aperture size of the Bragg cell, and the acoustic 
attenuation experienced by the acoustic wave as it transits 
the cell. This thesis 1S an attempt to provide a more 
detailed model by incorporating these factors ina first 
order simulation to aid in the study of acoustooptic spectrum 


and excision techniques. 
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II. ACOUSTOOPTICAL SPECTRAL ANALYSIS MODELING 


mee ACOUSTOOPTICAL INTERACTION THEORY 

As an acoustic wave propagates in an isotropic material, 
it produces periodic modulation of the index of refraction 
by means of the elasto~acoustic effect. The interaction of 
incident light waves with acoustic waves is termed acousto- 
optical coupling and is traditionally modeled utilizing 
coupled wave equations. The basic theory of acoustooptical 
interaction in isotropic media was well developed prior to 
the advent of the laser. An excellent discussion of early 
theoretical work can be found in Ref. 6 by Born and Wolf. 
The advent of the laser stimulated research in acoustooptic 
interactions as devices were developed for the deflection 
and modulation of incident laser light. Development of 
Superior acoustooptic materials and highly efficient broad- 
band transducers have given impetus to rapid advances in 
these acoustooptic applications. 

Various models describing acoustooptic interaction have 
been proposed. A thin grating model (as it assumes a negli- 
gible depth for the sound field) was proposed by Rosenthal 
[Ref. 7], Korpel and associates [Ref. 8], and Lambert [Ref. 9]. 
The thin grating model is essentially a first order solution 
of the diffracted light problem. A second approach known as 


the physical optics model incorporates the Green's function 
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integral, as proposed by Gordon [Ref. 10]. An extremely 
rigorous model, the mathematics involved are quite complex. 

A third approach is termed the phenomenological geometrical 
optics model as proposed by Johnson [Ref. 11] and ranks midway 
in complexity between the thin grating and physical optics 
models. 

For the purposes of this work a highly satisfactory 
approach to acoustooptical diffraction is found in a series 
of works by Chu, Kong, and Tamir [Refs. 12-16]. They have 
considered diffraction of arbitrary profiled optical beams 
by periodically modulated layers and have achieved closed 
form solutions for the resulting diffracted far-field 
intensity patterns. A summary of their results is presented 
below. 

As shown in Fig. 3, a bounded beam is incident ona 


periodically modulated dielectric layer with permittivity of 


the form 
€(z) = €2(1 + M cos 212/d) (1) 
where 
E2 = the relative permittivity of the slab in the 
absence of modulation 
M = the modulation index 
ad = periodicity 


The slab has thickness L and is bounded by a dielectric 


medium with relative permittivity ¢«, for xS0 and ¢3 for x-L. 
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The electric field of an incident beam may be represented as 


Enc (%, 2) = fo, exp ice. + Bo z) | a8 (2) 
with 

Bo = (27/rA) VE Sin 8 (3) 

e{l) _ (an/ryJer cos 8 

0 _ 1 Os 


yf (20/4)? Ein 309 (4) 


G(B.) is the angular spectral amplitude of the incident beam 


at the entrance plane x=0. 


oo 


a _- 
G(B8.) Lie | Pane (0,2) exp (-] B 2) dz 
=F1E,. (0,2)} (3) 
Letting 
eo (ORZiwa=obiz)sexp »(gbz) (6) 
with 
kg ; els) 
b = (27/A) VE. Sin 8 
= (20/4) Vex sin 0?) (7) 
where 
git) = angle of incidence of beam axis 
>? = angle of refracted beam 
F(z) = beam-profile at x=0 


then Equation (5) becomes 
G(Bo) = # [rw exp [-j (8 -b)z]dz (8) 
Pe: 





Various beam profiles F(z) and their corresponding angular 
spectral functions G(B.) are listed in Table I. 


For an optical beam at exactly Bragg incidence we have 


Ze = 0.) and 6{7) = of), ana 
be= (2nN/A)\/e, Sin - 
= (2n/A)Vex_ sin 6,°) 
= t/d (9) 


The far field pattern for the Bragg deflected beam 1S given 


by 
mo) = G(o) T(o) (10) 
where 
$ = (2m/r) Vez Sin 6 + t/a 
v(>) = 1s the complex transmission coefficient. 


A description of this coefficient may 
be found in Appendix C. 


Taus the far-field pattern for a Bragg deflected beam is 
essentially the Fourier transform of the anerture field at 


the boundary of the layer modified slightly by a transmission 


actor. 

While the results obtained by Chu, Kong, and Tamir are 
extremely useful in approaching the construction of a model 
for the direct detection acoustooptical spectrum analyzer, 
they are predicated on the periodic modulation of the layer 


by a Single frequency acoustic wave. 
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PeeeACOUSTOOPTIC SPECTRUM ANALYSIS THEORY 

The acoustooptic spectrum analyzer consists of only four 
components; a light source, Bragg cell, transform lens, and 
a photodetector array. The basis for acoustooptic spectrum 
analysis is the frequency dispersion of acoustooptic diffrac- 
tion. Essentially each frequency component of the input 
Signal will cause a deflection of the incident light beam. 

A typical interaction geometry is shown in Fig. 4 [Ref. 17]. 
A summary of the description of this interaction geometry 
from Ref. 17 is presented. 

Two discrete frequencies are shown with corresponding 
acoustic wavelengths A; and Az. These propagate in the x 
direction with velocity Wee The medium is optically iso- 
tropic with an index of refraction n. The thickness of the 
ultrasonic grating along the x axis is L. A collimated 
beam of light of freespace wavelength A is incident on the 
sound field at an angle os from the z axis in the xy plane. 

A sound wave at frequency f interacting with the incident 
light wave will generate a principal diffracted beam separated 


from the incident beam by an angle Os + 8. by the well known 


d 
Bragg law: 


_ dames Nite 
= (#-) (11) 
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where oO. + B 5 <- *0:. ) aeaid 


External to the medium 


ee i Se aes 
Oo. ~ 8 2S (a7) 


~ eee 
a ta (2) 


Thus the diffraction angle is approximately proportional to 
the input frequency. Figure 5 shows a classical optical ray 
trace diagram for the acoustooptical spectrum analyzer con- 
figuration. Figure 6 shows the corresponding optical plane 
schematic. In the output transform plane, position is 
approximately related to frequency by 


= =(z J (13) 


Ss 


where 

F = transform lens focal length. 

At this point, by considering the input signal as a 
composite of frequencies each of which will cause a deflected 
beam in accordance with the presented interaction theory, 
an approach to the construction of a model may be undertaken. 

Schulman [Ref. 18] and Vatz [Ref. 19] have proposed a 
set of optical block diagram symbols which show the process 
elements and indicate the mathematical function each element 
performs. Figure 7 shows the AOSA using these process 
elements. 


The laser light source §S and beam expander lenses, T 


ad 


and Toe are shown as providing an input to the Bragg cell 


M 
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Figure 5. Classical Ray Trace Diagram for an Acoustooptical 
Spectrum Analyzer 
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Figure 6. Optical Plane Schematic for an Acoustooptical 
Spectrum Analyzer 
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with transverse light distribution a(x) of general form, an 
arbitrary profile modified by spatial noise, lens aberrations, 
and Bragg cell aperture. The input signal f(t) 1s converted 
to an acoustic wave in the E/O block which represents the 
piezo-electric transducer. This acoustic traveling wave 
f(x-vt) propagates through the interaction region. 

In the Bragg cell interaction the light modulation can 
be considered a multiplication process in which the diffracted 
light power iS proportional to the product of the light inten- 


sity and the acoustic wave power. 


i(x) = a(x) £(x-vt) (14) 
where 
x = vee (253 


The transform lens T produces the Fourier transform. 


(Ss) 


3¥f{a(x) £(x-vt) } 


Il 


A(s) * F(s) exp [j2T7 S(v.t) ] (16) 
The transform is the convolution of the Fourier transform 

of the illumination function and the transform of the acoustic 

wave. The exponential factor 1s a Doppler frequency shift. 

As the detector 1S a Square law device, the output from the 


AOSA is the square of the magnitude of the intensity function. 


|A(s) * F(s)|- Cai) 

It has been shown that the output of the AOSA is the 
Square of the magnitude of the intensity function. The in- 
tensity function may be considered as the convolution of the 


Spatially transformed illumination function and the spatially 


Sa. 





transformed acoustic wave. Hecht [Ref. 20] treats the illumi- 
nation function as a spatial weighting function composed of 

a number of factors including optical amplitude profile, 
aperture size, and acoustic attenuation. The AOSA frequency 
plane intensity distribution would thus be proportional to 

the squared amplitude of the Fourier transform of the weighting 
function with respect to the spatial frequency. 


F(s) = f W(u) exp (-j2msu)du (18) 


2 if 
(ff - 3] 


finite period of signal analyzed 


where u 


x 


ft = the normalized frequency 

Further refinement in the basic model is made by the 
consideration of the effects of transducer efficiency and 
diffraction efficiency. These factors modify the frequency 


components in the acoustic wave. 


C. AOSA MODEL CONSTRUCTION 


1. Irradiance Profiles 
The choice of the optical beam profile incident on 
the Bragg deflector is a complex decision. It would appear 
desirable to fill the entire deflector aperture with a uniform 
intensity light below the saturation level of the AOSA. The 
irradiance profile in this case would be that of a plane wave. 
A uniform intensity filling the entire aperture will 


maximize acoustooptic interaction and provides maximum 
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incident light to be contributed to the deflected beam 
components. Abrupt windowing will truncate the beam profile 
and cause sidelobes to appear. These sidelobes, if they are 
of sufficient magnitude, may cause erroneous components to 
appear at the image plane or mask low intensity components. 
To reduce the sidelobe levels, a less discontinuous 
transition in the absorption-transmission regions is necessary. 
A suitable choice would be a GausSian irradiance profile with 
a spherical phase wavefront. This irradiance profile common- 
ly occurs as a laser output and has well understood propaga- 
tion behavior. Smith includes a handy propagation program as 
part of his work [Ref. 3]. This program calculates Gaussian 
beam phase front radius and spot size for propagation through 
any number of thin lenses spaced at arbitrary distances. 
Since the Gaussian incident beam profile would not fill the 
entire aperture with uniform intensity and assuming continued 
operation below the AOSA saturation level, less light inten- 
Sity Would be available for diffraction into representative 
beams. This will result in the loss of sensitivity and the 
possibility of an inability to detect low amplitude frequency 
components in the presence of strong frequency components. 
Other beam irradiance profiles are possible and 
might be employed for their specific deflected beam charac- 
teristics. Figure 8 taken from Ref. 15 shows the deflected 
beam profiles for various incident profiles. Such curves 


are useful for comparison of expected sidelobe levels. 
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Hecht's [Ref. 5] implementation of the Gaussian 


irradiance profile is of the form: 


wi(t,t,T) = exp aaa - 5) | (19) 
where 

T = D/2W 

D = effective aperture 

yy = Gaussian waist size 


A significant design question in the utilization of Gaussian 
irradiance profiles is the choice of the ratio of the aper- 
ture to the ae intensity. This is discussed also by Smith 
[Ref. 3]. Beam spot size is defined as the perpendicular 
distance from the beam centerline to the point where the 
irradiance has dropped by a factor of Page To use the full 
information capacity of the Bragg cell, the spot size must 
be enlarged to illuminate the complete cell aperture. [In 
this situation the aperture will severly truncate the Gaussian 
beam resulting in a profile similar to that of a sinc 
function. Varying the ratio of the spot size to the aperture 
will alter the level of the sidelobes of the profile. Smith 
limited his spot size to avoid sidelobes and thus isolate 
their effect in his study of optical excision. 

Hecht expresses his spatial weighting functions in 
equivalent time variables. The linear relationship between 
the spatial and time variables and the analogous relationship 


between the spatial frequency and frequency provide great 
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Mmeexebility in the construction of the computer program to 
model the AOSA. The Fourier and inverse Fourier transforms 
utilized in the AOSA program make no distinction as to the 
transforming variables. 

Only plane and Gaussian irradiance profiles were 
included in the AOSA program. It would be a simple matter 
to incorporate other profiles such as the Lorentzian or 
triangle beams. However, the frequency with which the more 
exotic profiles would be encountered does not merit the 
addition. Figure 9 shows the comparative image plane 
profiles for both irradiance profiles with a monotonic 70 MHz 
input Signal. In Figure 9 the plane wave illumination chosen 
completely illuminated the aperture of the Bragg cell. The 
GausSian illumination used an arbitrarily chosen 3.8 milli- 
meter spot size. Comparison shows a near 50% reduction in 
the level of the Gaussian sidelobes to the level of the plane 
wave sidelobes but with a tradeoff of a decrease in the peak 
intensity response. The units for the Image Plane Intensity 
axis in this figure and all succeeding AOSA program figures 
are arbitrary. 

Figure 10 shows the image plane intensity for an 
incident Gaussian beam with a spot size of 1.945 millimeters. 
This is the “optimum" spot size according to Smith's 
criteria for sidelobe suppression. The sidelobes are barely 
visible but again a tradeoff results in an additional decrease 


in peak intensity. Use of a Guassian beam with this spot 
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size narrows the portion of the Bragg cell aperture receiving 
full illumination. The portion receiving the full illumina- 
tion is termed the effective aperture. Its effect will be 
discussed in a later section. Figure 1l shows for comparison 
the image plane intensity profile of a fully illuminating 
plane wave. 
eee aperture Size 

The universal feature of realistic weighting functions 
is that of a finite aperture corresponding to a finite period 
of signal t represented by the unity amplitude window function 
[Ref. 5]. 

w(t,t) = rect (: - :) (20) 

18 2 

where 


T= D/V . 


Early approaches to acoustooptical spectrum analysis modeling 
considered the input signal aS a continuous wave. The effec- 
tive aperture of the selected Bragg deflector truncated the 
incident light, such that only the segment of acoustic wave 
within the optical aperture region was intercepted. That 
segment of acoustic wave corresponds to the finite period of 
Signal being analyzed. 

Improvements in Bragg deflector manufacturing tech- 
niques have yielded high quality, low acoustic velocity 
deflector crystals. These crystals are of increasingly larger 
dimensions. The Rayleigh limit resolution, N, is defined as 


the total deflection angle divided by the angular spot size. 


39 





_— 
— 
Ch) 
= 0.20 
Be! 
— 
= 0.15 
Le 
z= 
= 0.10 
| 
oa 
= 60.05 
G) 
(eu 
> 
= 0.00 
Eagure 11. 


sl), =) Ue Jide elicle 90. 


eee) tia) 


Bales) | BINGE CO) 


Image Plane Intensity Profile (Magnitude) Produced 
by a Monotonic 70 MHz Signal with an Incident Plane 
Wave and a Simple Transducer 


40 





The relation simplifies to 
N wae CC) oe Ae (2) 


where 


eA) acoustic transit time across the optical aperture 


bids 


deflector bandwidth 
For the AOD-70 deflector (assuming that full aperture illumi- 
nation is provided) 

t, = 2.56 x D (microsec) (22) 
where D is the optical aperture in centimeters. 

Using the AOD-70 parameters of a 40 MHz bandwidth and a maxi- 
mum t, of 10 microsec based on the full available aperture 

of 3.9 cm, calculation yields a Rayleigh resolution capacity 
of 400 spots. To fill the full aperture anamorphic optics 
such as cylindrical lenses or telescopes must be used. Using 
standard procedure for the detector array utilization each 
spot 1s focused to fall across two adjacent photosites. The 
Rayleigh resolution capacity of 400 spots provides for a 
maximum frequency resolution of .1 MHz for the full 40 MHz 
bandwidth. 

Using a Gaussian beam with spot size of 1.945 milli- 
meters to fully suppress the sidelobes in the image plane 
intensity profile results in an additional penalty ina 
decrease in the Rayleigh resolution capacity of the "time- 
bandwidth product," as given by Equation 21. The interaction 
length is no longer the Bragg cell aperture but rather the 


effective aperture, in this case estimated as the region 
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twice the spot size symmetrical about the cell spatial midpoint. 
This results in a smaller optical aperture, a smaller acoustic 
transit time, and a resultant Peeeeaee in the Rayleigh resolu- 
tion. Thus, while a choice of spot size of this magnitude 

will suppress sidelobes, it requires a possibly unacceptable 
tradeoff in frequency resolution. 

The larger apertures bring an additional consideration 
to the AOSA model construction. A typical application input 
Signal pulse might be of one microsecond in duration while 
the total acoustic transit time for the AOD-70 is ten micro- 
seconds. A one microsecond pulse input to the deflector will 
persist in terms of the frequency spectra it contains and 
therefore, in the deflected beam components incident on the 
image plane for the full ten microsecond travel time. This 
effect of "time spreading" of the r-f signal is significant 
for its enhancement of the detection probability for short 
duration signals. In the application of AOSA techniques to 
Electronic Support Measures (ESM) receivers, this is of prime 
import. A 0.1 microsecond r-f signal, which is currently 
beyond the capability of many traditional technology base 
receivers, would be stretched to a 1.0 microsecond pulse 
greatly improving its detection probability. It should be 
noted that while this "time spreading" increases detection 
probability, it reduces the capability for accurate direct 
pulse width measurement, also a highly desirable factor in 


ESM application. 
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In summary the aperture size weighting function 
imposes three constraints on the AOSA model. The physical 
parameter of the optical aperture is modeled by a simple 
unity amplitude window function. The Rayleigh resolution 
and thus frequency resolution is simulated by the technique 
of zero padding while applying the Fourier transforms to 
achieve a digital record length which will yield the appro- 
priate frequency resoltuion. Zero padding is a digital 
implementation technique used in Fourier analysis whereby 
the desired sampled signal data is entered followed by a 
data string of zeroes. This increases the record length to 
provide a desired resolution in the transformed domain with- 
out affecting the sampled data validity. The transformed 
domain resolution is simply the inverse of the record length. 
The finite signal length is modeled with a unity amplitude 
window function implicit in the input signal as it is digitally 
implemented. Since the period of signal analyzed may be 
shorter than the acoustic travel time of the optical aperture, 
1t 1s desirable to analyze the deflection caused by inter- 
action of the pulse acoustic segment at various points as it 
traverses the optical aperture. A shiftable digital array 
was utilized controlled by a time delay variable; zero delay 
corresponding to a symmetric straddling of the acoustic pulse 
Segment about the spatial midpoint of the optical aperture. 
This time delay variable forms the basis for a dynamic AOSA 


program to be discussed in a later section. 
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3. Acoustic Attenuation 

As the acoustic wave travels the length of the Bragg 
deflector, it is subject to acoustic attenuation. The ampli- 
tude of the acoustic wave decreases, with a resultant decrease 
in the ability to interact with the incident light wave. 
Hecht [Ref. 5] treats the acoustic attenuation as a spatial 
weighting function of the form; 

w(t,a,tT) = exp(-at) (23) 
where a is a frequency dependent loss factor 

ages 


Qq = (24) 
20 109) e 


with 


: dB 
re = usec GH™ 


Values for oe have been tabulated for common acoustooptic 
modes. Acoustic attenuation typically follows power law 
frequency dependence with exponent, m, near 2. 

The AOD-70 Bragg deflector is manufactured with high 
quality flint glass as the interaction medium and Indium 
bonded Lithium Niobate for the piezo-electric transducers. 
The Oo values for this type of deflector range from .1 to .15 
(dB/microsecond GHZ*) . A plot of the more stringent condi- 
tion is shown as Figure 12. Calculation of the acoustic 


attenuation for the AOD-70 device parameters in terms of 


a normalized transmission factor was carried out. The 
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results are plotted as Figure 13. These results show a near 
unity transmission for the region of concern, from 50 to 
90 MHz. Since the transmission loss was negligible for this 
application, an acoustic attenuation factor of Hecht's form 
was not implemented in the AOSA programs. 

If the use of the AOSA programs iS envisioned in 
higher frequency ranges, acoustic attenuation will become a 
Significant factor. Figure 13 shows increased loss as the 
frequency nears 1 GHZ. Further discussion of the effect of 
acoustic attenuation on weak truncations (Gaussian beam) and 
strong truncations (uniform optical beam) is given by Hecht 
fRef. 21]. 

Paolo heaction ELLiciency 

For a simple transducer the Bragg angle cannot be 
exactly matched by the incident collimated light beam for 
all frequencies over the bandwidth of concern. As a result, 
for Bragg regime deflection the intensity of the diffracted 
beam is a function of the applied frequency. The ratio of 
the diffracted beam intensity to the intensity of the incident 
beam is defined as the diffraction efficiency. An analysis 
of Bragg regime deflection was conducted by Reagan in Ref. 2 
and is summarized here. 

Hecht [Ref. 5] expresses diffraction efficiency in 


a simplified form as: 


T 
1 : 25 oe (FFF?) | 
as See (72) SA eS = (25) 
O E r (FEAF ) 
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with 


n= Taft va” (26) 
22 ov® 4 cos’ 6 

and 

= ale etek ae (27) 
where: I, = diffracted beam intensity 

I, = undeflected beam intensity 

ve = free space laser wavelength 

n = refractive index 

p = elastooptic coefficient 

O09 = mass density 

An = acoustic wavelength 

Vv = acoustic wave velocity 

E = acoustic power 

H = acoustic beam height 

6 = Bragg deflection angle 

L = interaction length 

La = characteristic interaction length 

F = normalized center frequency = l 

a = input frequency value normalized to the center 


frequency value fy 
Equation 26 shows the quantity n to be constant with 
respect to the input signal frequency. Thus the initial sin 
Squared term of Equation 25 may be ignored when considering 
the bandshape of a normal acoustooptic beam deflector. A plot 


Of Equation 25 for various L/L, ratios is shown as Figure 14. 
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The 3dB bandwidth is given by Equation 28 and is also from 
Pet. 5. 


ho 1a fy L/L (28) 


For application of a Bragg regime deflector as an 
acoustooptical spectrum analyzer, it is desirable to maintain 
a bandshape aS near unity and as uniform as possible over 
the bandwidth of concern. The characteristic interaction 
length is shown by Equation 27 to be constant with respect to 
dimensions of the Bragg deflector. The trend for AOSA 
applied Bragg deflectors is that of increased size with an 
attendant increaSed interaction length L. This reflects 
both improvement in manufacturing techniques permitting the 
production of high quality large deflectors and the desire 
for increased resolution made possible with larger inter- 
action length. The increased interaction length, however, 
results in higher value L/L, ratios with the observed 
effect of narrowed bandwidths. 

Attempting to increase bandwidth by shortening the 
interaction length is of limited benefit since the diffraction 
shifts from the Bragg regime to the Raman-Nath regime. 
Raman-Nath regime diffraction is unsuitable for this appli- 
Cation as it supports many higher order diffraction modes 
resulting in interference created by additional beams. 
Increasing the refractive index or elastooptic coefficient 
of the interaction medium, decreasing the acoustic beam 


height, or decreasing the acoustic velocity or the mass 
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density of medium can improve the diffraction efficiency to 
a limited extent. An increase of acoustic beam power will 
increase diffraction efficiency. The increased power also 
tends to heat the interaction medium, creating temperature 
gradients which degrade resolution and information capacity. 
Increased difficulty in internal reflection damping is 
experienced, permitting unwanted Signal persistence. 

The only effective manner in which to increase the 
bandwidth is to increase the cell's center frequency f5° 
The maximum frequency which the input transducer can effec- 
tively couple to the medium is determined primarily by the 
bonding technique used to attach the piezo-electric strips 
to the interaction medium. Thus the ability to increase the 
center frequency is dependent on the input Signal transducer. 
Normally Bragg deflectors are limited to operation over one 
Sccave . 

In summary, for a simple transducer little can be 
done to improve diffraction efficiency or increase bandwidth. 
Beam steering, however, can significantly improve both 
diffraction efficiency and bandwidth. Commonly employed 
beam steering techniques include first order planar and 
stepped arrays. Figure 15 shows a typical first order 
Stepped array. Pinnow [Ref. 22] discussed the relative 
merits of simple, planar and stepped array transducers. 

The larger dimensioned Bragg cells deflectors are almost 


always of the stepped array type. The AOD-70 deflector is 
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a five element stepped array. Figure 16 is the manufacturer's 
supplied efficiency versus frequency response curve. 

The AOSA program has both a simple transducer and 
a stepped array approximation option. The simple transducer 
response is computed in accordance with isotropic theory 
based on the AOD-70 device full dimensions. This permits a 
comparison with the actual stepped array response. Figure 17 
shows the computed diffraction efficiency curves for various 
L/L, ratios. These were obtained by sweeping a single varia- 
ble frequency over the band of concern and extracting the 
magnitude value for the instantaneous frequency applied. In 
comparison with Figure 14, which contains the bandshapes for 
the diffraction efficiency theory, a rippling is noticeable 
at the lower frequencies. This rippling is an undesirable 
@rtitact of the digital Fourier transform application. The 
ripple amplitude is comparatively small and falls outside of 
region of concern; hence it was not corrected. 

The stepped array bandshape that was included in the 
model iS an approximation based on the near uniform response 
that was experimentally measured for the AOD-70 device. A 
Slight roll-off in response occurs at the endpoints of the 
band, but for a first order approximation a uniform window 
function utilizing the mean efficiency was employed. 

To compare the response of the models of the two 
transducers, an input signal consisting of an in-phase 


frequency comb with 5 MHz spacing was applied. Figure 18 shows 
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the simple transducer response. Figure 19 shows the stepped 
array response. Comparison shows severe attenuation of the 
frequency components beyond 80 MHz and a lesser degree of 
attenuation from 50 to 60 MHz for the simple transducer. 
On the other hand the stepped array transducer shows a very 
uniform response but with a tradeoff in decreased overall 
intensity magnitude. 
oe Leanisducer Efficiency 

The measure of the transducer's ability to convert 
the r-f input Signal components to representative acoustic 
waves 1s termed the transducer efficiency. The transducer 
efficiency is dependent upon the transducer material, the 
bonding technique used to fasten the transducer to the inter- 
action medium, and the transducer impedance matching. A funda- 
mental tradeoff relationship exists between the flatness of 
the frequency response and the amount of power delivered by 
the transducer within the desired band. The power levels 
used to drive the transducer are relatively low, thus we can 
afford a power conversion inefficiency and seek as uniform 
a response as possible. All transducers operate at less than 
unity efficiency; however, for purposes of the model, it was 
assumed that the desired uniform response was of paramount 
importance. An ideal case would provide a transducer efficiency 
response with little power attenuation and a uniform amplitude 
for all frequencies within the range of the system application. 


Current technology has achieved near uniform response bandshapes 
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over bandwidths of 70 MHz, as is present in the AOD-70 cell, 
Chang [Ref. 23] discusses the state of current transducer 
technology. 

Based on the above the AOSA model was implemented 
with a unity amplitude window for the transducer efficiency 
function. The program structure permits easy input of non- 
unity bandshapes should this be desired. Extreme broadband 
applications of acoustooptical spectrum analysis will result 
in nonuniform frequency responses requiring such modifications 


before using the AOSA program. 


D. AOSA MODEL UTILIZATION AND EXTENSIONS 

The AOSA program was designed to model direct detection 
acoustooptical spectrum analyzer systems in general. [In its 
current configuration it utilizes the component parameters 
of the NPS laboratory system with the addition of the AOD-70 
Bragg cell. The coding for the Fortran program was selected 
to maintain as much generality as possible, so as to permit 
application for analysis under other choices of components. 
An attempt was made to minimize the number of preliminary 
calculations a prospective user would have to perform prior 
to utilization of the program. For example, the program 
computes the Fourier transform lens focal distance given the 
deflector physical parameters, the laser characteristics, 
and the desired detector array spot size. A list of the 


required input physical parameters 1s given in Table II. 
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TABLE II 


AOSA Modeling Parameters 


DCELL = effective optical aperture of Bragg cell 

LW = Bragg cell interaction length 

vo = acoustic velocity of interaction medium 

REFIND = refractive index of interaction medium 

LAMBDA = laser wavelength 

BW = beam width of collimated Gaussian beam (3 dB) at 
the Bragg cell aperture 

SPOT = desired detector array spot size 

PULSE = r-f signal pulse duration 

TDELAY = time delay for examining signal pulse at varied 


points during the acoustic transit period 
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The remaining program variables are primarily concerned with 
control of the program including output formats and the var- 
ious subroutines utilized within the program. 

The main exception to the above is the created variable 
TDELAY, for time delay, which permits the user to examine 
the representative acoustic pulse at frozen intervals of time 
during its transit of the optical interaction aperture. MThis 
Suggested the possibility of adding another dimension to the 
AOSA program, to create in effect a dynamic AOSA model which 
would examine the response of the acoustooptic spectrum 
analyzer periodically by looping through the basic static 
algorithm with decreasing time delays. The results are 
presented in a three~dimensional isometric or perspective 
projection. This permits the user to simulate the response 
of the image plane intensity as the acoustic pulse transits 
segments of the optical interaction region rather than at a 
Single fixed point of time. To simplify the graphic output 
the present configuration of the dynamic AOSA program (called 
TAOSA) is set to examine the first fraction of a microsecond 
as the acoustic pulse enters the optical region. Computer 
Simulations were conducted over the entire transit period 
and to verify the symmetry of response as the acoustic pulse 
exited the optical region. 

Both the static and dynamic programs provide the user 
with options to simulate various system configurations. 


Types of transducers, incident beam profiles, and output in 
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either relative magnitude or decibel form is available. The 
static program offers a choice of several levels of graphical 
output, from a coarse line plot to two Versatec offline 
plotting options. Currently both programs utilize a deter- 
mMinistic form for the input r-f signal. The programs could 
be easily modified to accept sampled r-f data. 

The following series of figures are typical of the outputs 
of the AOSA and TAOSA programs. They were selected to demon- 
Strate the ability of the programs to model a variety of 
complexity in input Signals, to demonstrate the quality of 
resolution achievable by an acoustooptical spectrum analyzer 
and for comparison with the known spectral content of the 
inputed signals. 

Figure 20 is the input r-f signal magnitude for a 70 MHz 
monotone. Figures 21 and 22 are outputs using the decibel 
option of the AOSA program. Figure 21 shows the relative 
Sidelobe levels by interaction of the monotone with an inci- 
dent plane wave. For comparison Figure 22 shows the decrease 
in the relative sidelobe level produced by the monotone with 
an incident Gaussian beam of 1.045 millimeter spot size. 
Figures 23 and 24 show the dynamic response of the same 
monotone with the two wave types. The x axis on all TAOSA 
program graphics represents time as the acoustic pulse transits 
the Bragg cell aperture region. Time zero corresponds to the 
instant when the leading edge of the acoustic wavetrain 


crosses into the aperture region. Figure 23 shows the dynamic 
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Figure 20. Input RF Signal Magnitude for a 70 MHz Monotone 


63 





= —e 


lt 2 


—  —————_— 


(DB) 


_— 

— 

m= 20. 

Y) 

= 

Li] 

— 

= -ud. 

= 

Li 

Zz 

ae 

m6 - 0. 

a 

EL 

G) 

= -3gd. 

> 

— 
Figure 21. 





be 


, 


wall | 


Piece.) (Mid 7) 


=) OF 


BLS aliG ee) 


Image Plane Intensity Profile (Decibel) Produced 
by a Monotonic 70 MHz Signal with an Incident 
GausSian Beam and a Stepped Array Transducer 


64 | 





(DB) 


_— 

— 

= 6 C- 20. 

ws) 

rae 

Li} 

— 

= -40. 

Ey 

ai 

a -cO. 

Qu ; : 

= | if 

OC | 
au -90. ia! Li. 

= SO. 60. 70. 80. a0. 


eect) | ih 2) 


SISTANCE . COs 


Figure 22. Image Plane Intensity Profile (Decibel) Produced 
by a Monotonic 70 MHz Signal with an Incident 
Plane Wave and a Stepped Array Transducer 


65 





Ta 
a) 
‘ 

a 





ale 
ww 
¢ 


Law44 usec 


NOS AN : 
S 
‘ AD s 
¢ : ‘ : 
a0 FANN 
ah 
BD sae Hes 443 
‘ ANAS 
* * A 
YAS AGS SOUR NOAA ae 
< we a Y 
He “Ne 
e' es &* « Lan 
8 - fl 
As © * ? 
et Ne A Gi 3 
. q ee . We 2 4 
’ SoA 


* 5 
Z < 

We'd 
‘Q 


AWN 





UNITS 


(usec) 


Time 


X axils 


: 

Oy 
OHM 
OU a 
rages) 1S 
Ovy 
Beal 
ts! ceo 
0 DH 
“Gd 
FP) 
11 
nw 
“4 4 
dee 
Oo 
SON 


Produced 


(Magnitude) 
by a Monotonic 70 MHz Signal with an Incident Plane 
(Dynamic Model) 


Wave and a Stepped Array Transducer 


Image Plane Intensity Profile 


megure 23. 


66 





response of the monotone with an incident plane wave. 
A much slower response is seen in Figure 24 showing the 
response produced by an incident Gaussian wave of 1.945 
millimeter spot size. The Gaussian response is shown to 
occur later in the cell aperture transit period and to last 
for a relatively short time when compared to the plane wave 
response. The symmetric roll-off is already visible at the 
rightmost edge of the figure. 

Figure 25 is the input r-f signal magnitude for a simple 
Sinusoidally modulated FM signal with 8 = 8.0 and fa = 5 MHz. 


8 is the modulation index and is defined as 
ate (29) 


Figure 26 shows the image plane intensity profile produced 
by the interaction of the signal shown in Figure 25 with an 
incident plane wave. Figure 27 shows the image plane inten- 
Sity profile produced by the same signal interacting with an 
incident Gaussian beam of 1.945 millimeter spot size. 
Figures 28 and 29 are the analagous dynamic program outputs. 
Again the Gaussian beam interaction shows a slower and 
smaller response. The advantage of sidelobe suppression, 
however, iS apparent in the extremely clear presentation of 
the FM intermodulation products. 

Figure 30 shows the image plane intensity profile produced 
by a linear chirp with an incident plane wave. Figure 3l 


shows the profile produced by the interaction with a Gaussian 
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beam. For both the chirp parameters are ‘ff = 80 MHz, 

fy = 50 MHz, and t = 2 microseconds. Both show hich and 

low frequency ripples. The high frequency ripple is an arti- 
fact of the digital implementation of the linear chirp. The 
lower frequency ripple present is the effect of the funda- 
mental and sidelobe interaction. The Gaussian interaction 
shows the effects of its inherent slower response and degraded 
frequency resolution in the attenuation of magnitude over the 
lower halfband frequency components. 

Figure 32 shows the dynamic response of the same linear 
chirp and an incident plane wave. The severity of the inter- 
action of the fundamental and sidelobe components, as addi- 
tional frequency components enter the Bragg cell, manifests 
itself in a dynamic transient response. Only after a con- 
Siderable interval does the response settle into a steady 
state displaying the expected chirp spectral profile. For 
comparison Figure 33 shows the same linear chirp interacting 
with an incident Gaussian beam of 1.945 millimeter spot size. 
In this case the sidelobe suppression is a distinct advantage 
producing a smooth steady state response almost from the 
onset of the interaction. It must be kept in mind that the 
interaction will continue to begin much later than that of 
the plane wave and thus the plane wave will achieve steady 
state earlier than the Gaussian. 

The TAOSA program also has the decibel output option. 
Figure 34 shows the relative decibel image plane intensity 


of a monotone of 70 MHz interacting with a plane wave. 
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III. SUMMARY AND CONCLUSIONS 


me. DUMMARY 

This thesis has presented an analysis of Bragg diffrac- 
tion theory and direct detection acoustooptical spectrum 
analysis modeling techniques. Based on the study of the 
various modeling theories a suitable approach for a first 
order simulation was selected and used to develop a Fortran 
program modeling the direct detection acoustooptical spectrum 
analysis process. Extensions of the basic theory have been 
added to include the factors of diffraction efficiency, 
transducer efficiency, irradiance profiles of the incident 
laser ination , the aperture size of the Bragg cell, 
and acoustic attenuation of the generated acoustic wave. 
A second program dynamically models the AOSA process incor- 
porating a time variable to simulate the response of the AOSA 
as the acoustic wave transits the interaction region. The 
programs were used to model the NPS laboratory AOSA using 


the parameters from that system. 


B. AREAS FOR FURTHER STUDY 
Further effort related to this thesis is required in the 
following areas. 
1. Nonlinearities 
The theoretical basis for the AOSA model program does 


not take into account the nonlinear effects of crossmodulation 


80 





(the compression of one signal due to another) or inter- 
modulation (the interaction of two or more signals). While 
the restriction of the operation of the AOSA to the Bragg 
regime eliminates much of the effects of these nonlineari- 
ties, the overall effect on the replication efficiency is 
still at issue. Further analysis is needed to determine 
if these factors seriously degrade the model's accuracy. 
Hecht discusses this issue in Ref. 21. 
2. Time Varying Signals 

The ability of the AOSA to respond to rapid time 
varying signals (which will depend on the speed of response 
of the detector array and the Bragg cell) needs to be 
examined. Currently the TAOSA program provides for a fixed 
delay for the simulation of overall system response time. 

3. System Components 

The physical response and performance of the AOSA 
systems components warrant further study. Currently only 
a gross Simulation of the system response time can be made 
by insertion of a fixed time delay. Closer examination of 
the speed of the detector array with particular emphasis on 
the uniformity of photosite response is needed. Integrated 
optical devices for AOSA and AOSE applications will use 
semi-conductor lasers. These lasers will need to be 
examined for the system effects produced by fluctuations in 
intensity and variations in coherence. The current HeNe 
laser used in the NPS system should be examined for its 
intensity and coherence fluctuations and resultant effects. 
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c. CONCLUSIONS 

Comparison of the AOSA model image plane intensity 
profiles with known AOSA system output spectral profiles 
demonstrates the AOSA program as a satisfactory first order 
Simulation. The acoustooptical spectral analysis process 
yields high quality replication of the spectral content on 
input signals. Future coupling of the AOSA program with 
optical excision routines previously developed may provide 
additional insight into the application of acoustooptical 
Spectral analysis and excision for use in wideband, high 


density signal environments. 
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APPENDIX A 


ACOUSTOOPTIC SPECTRUM ANALYZER (AOSA) PROGRAM 


This Fortran program calculates the relative Bragg 
deflected image plane intensity profiles for multifrequency 
input r-f signals to a modeled acoustooptic spectrum analyzer 
(AOSA). The physical parameters used in the modeling of the 
AOSA are discussed within the program and are summarized in 
Table II of the thesis. 

The program provides several additional selectable 
variables: 

ITEST - provides for graphical and numerical or 
Graphical only Output 

IWAVE - provides for choice of plane or Gaussian 
incident light waves 

ITRANS - provides for choice of simple or first 
order stepped array transducer simulations 

IDB - provides for choice of graphical output in 


magnitude or relative dB formats 


IGRAF - provides for choice of line printer or 
Versatec graphical plotting 
IPLOT - provides for choice of PLOTG or DRAWP 


Versatec plotting subroutines 


Use of these variables is discussed within the program. 


8 3 





Additional alterable variables exist within the program 
related to the Fourier transform process. It is believed 
that the present configuration will suffice for most general 
applications. 

This program utilizes the PLOTG and DRAWP Versatec 
software subroutines, the line plot subroutine PLOTP, and 
FOURT, which calculates the forward and inverse Cooley-Tukey 
Fast Fourier Transform of multidimensional complex data. 
Prospective utilizers at the Naval Postgraduate School should 
refer to the W. R. Church Computer Center for detailed 


information. 
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AEEENDIX B 


DYNAMIC ACOUSTOOPTIC SPECTRUM ANALYZER PROGRAM 


The dynamic acoustooptic spectrum analyzer program (TAOSA) 
is essentially a variation of the AOSA Fortran program. It 
computes the time delay to shift the acoustic pulse outside 
of the optical interaction aperture and then incrementally 
decreases this delay to simulate the action of the pulse as 
it transits the optical interaction region. The system 
response time may be simulated by adding a fixed time delay. 

The variables used in the program are discussed both 
within the program and in Appendix A. An additional variable, 
TINT, 1s added as an acceleration factor that varies the in- 
crement size based on the level of intensity of the incident 
laser beam. 

The results are presented as a three-dimensional isometric 
Or perspective plot. Perspective users at the Naval Post- 
graduate School should refer to the W. R. Church Computer 
Center for detailed information on the three-dimensional 


plotting routine PLT3D1 in use. 
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INITIALIZE D ARRAY VALUES AT ZERO 


INITIATE LOOP TO PROVIDE MULTIPLE TIME VARIANCE FOR DYNAMIC MODEL 
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APPENDIX C 


ARBITRARY PROFILE SOLUTION DETAILS 
{after Ref. 15) 


The complex transmission coefficient T(¢) is defined as 
= a y _ ap & 
T(o) 4k 76 (o, a JCAL, Ba? exp ( jk, L)/Det (Cl) 


where 


Det = (a oe - ee op - 8 Bae 
= a = (C2) 
A1a2(A_, 2 ye aa Ba? 
l 4 
oS faa, {B2, — 85 + ((62, ~ 83)? + (Mbgen)7\0/"), (C3) 
If ° . 
Qty = nie Oo! — 83 — [(G2, — B3)2 + (Ma3e0)?]1/2), (C4) 
nibs Ria (, RL 
sens (4) (x BS 
X [exp (~tk3,L) — RORY exp(tk3,L)|, (C5) 
Boe =k3 [i + (1 + 
Se Ne ee ke 
x fexp(—ik3,L) — R3PR33 exp(ik3.L)], (C6) 
kt. —ké 
ef = atl = om 
Re: be + Re, p,o2a,b, tj = 1,2,3, (C7) 
ki. = Rov €, cod — V(2/\}* 6; a Ri (eo) 
hi, = (je, — 82.) 2 (ey 
RO, = \[1 — (Y)May|kjeq — 33}!/2, (C10) 
3, = \(L — ('4)Mag)Rgeq — 33]/2, (cine) 
RG, = (R5e€3 — (33)1/2 (Ei2)) 
Ro, = (Ries — B2,)'/2 (E13) 


M = index of Modulation of the permittivity of the Bragg 
cell medium 


~ =relative permittivity of air =1 
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relative permittivity of the Bragg cell medium 
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output aug@te 


2m1/A 
(20/2) ve sin 8 (C14) 
3 - 2%/d (C15) 
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